TeV neutrinos from accreting X-ray pulsars 



On 
O 
O 
(N 



Wlodek Bednarelfl 

Department of Astrophysics, University of Lodz, 90-236 Lodz, ul. Pomorska 149/153, Poland 

(Dated: May 30, 2009) 

Pulsars inside binary systems can accrete matter that arrives up to the pulsar surface provided that 
its period is long enough. During the accretion process, matter has to be accelerated to the rotational 
velocity of the pulsar magnetosphere at the distance where the balance between the pressure of 
matter and the magnetic field is achieved. At this distance, a very turbulent and magnetized region 
is formed in which hadrons can be accelerated to relativistic energies. These hadrons can interact 
with the very strong radiation field coming from the hot polar cap on the neutron star surface 
created by the in-falling matter. We calculate the neutrino event rates produced in an km 2 detector 
that can be expected from accreting millisecond and classical X-ray pulsars at a typical distance 
within our own Galaxy. 

PACS numbers: 97.80Jp,97.60Gb,96.40Qr,98.60Ce,98.70Rz,98.20Hn 
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I. INTRODUCTION 

A large number of neutron stars (NS) has been dis- 
covered as close companions of low mass and high mass 
stars (so called low mass and high mass X-ray binaries, 
see [H, 0] ) • These binaries are characterized by a strong 
X-ray emission that is produced as a result of matter ac- 
cretion from the companion star onto the neutron star. 
In some cases, the matter arrives up to the NS surface 
at the magnetic polar cap region which creates a very 
hot spot. This class of X-ray binaries shows pulsations 
correlated with the rotational period of the NS (so called 
X-ray pulsars, see @). X-ray binaries may contain a long 
period NS as well as a millisecond NS. The long period 
X-ray pulsars are usually found in isolated binary sys- 
tems whereas the millisecond pulsars (MSPs) are mostly 
contained within globular clusters Q. According to the 
standard scenario, the MSPs are old pulsars spinned up 
by the angular momentum transfered to the NS from the 
accreting matter. Many of them have been captured by 
stars whose density inside the globular cluster is very 
large. On the other hand, long period X-ray pulsars are 
found inside binary systems in which NS is formed in a 
supernova explosion of a more massive companion. 

The matter that is falling onto the NS as a result of the 
overflow through the inner Lagrangian point usually con- 
tains a large angular momentum. This matter forms an 
accretion disk that can be disrupted at the inner radius 
due to the pressure of the magnetic field of the NS. In the 
case of accretion from the wind of the massive star, the 
process occurs more spherically due to the isotropization 
of the accretion flow by the shock in the massive star 
wind. In some binaries the process of accretion ca be ad- 
ditionally stimulated by the re-radiation of stellar surface 
of the companion star by the X-ray radiation from the 
NS. Accreting pulsars can suffer spin up and spin down 
periods depending on the conditions inside the binary 
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system and the parameters of the MSP (i.e. the mass loss 
rate of the companion star, distance between the stars, 
neutron star velocity, surface magnetic field strength and 
the MSP rotational period) . These basic parameters de- 
termine the specific accretion phase of the matter onto 
the NS, i.e. whether it is in the accretor, the propeller, 
or the ejector phase The ejector phase is characteris- 
tic for a fast rotating NS whose magnetized wind is able 
to prevent any accretion below the light cylinder radius. 
In the propeller phase, matter can enter the inner pulsar 
magnetosphere. However, it is blocked by the NS fast ro- 
tating magnetosphere at a certain distance from the NS 
surface due to the centrifugal force. Then, most of the 
matter is expelled from the vicinity of the pulsar and the 
accretion onto the neutron star surface may occur only 
episodically when the pressure of the accumulated matter 
overcomes the pressure of the rotating magnetosphere. 

In this paper we are interested in the acceleration of 
hadrons during the accretion of matter from the compan- 
ion star onto the magnetized neutron star. Hadrons inter- 
act with the thermal radiation from the NS surface and 
produce neutrinos which might be observable by the large 
scale neutrino telescopes (IceCube, KM3NET). Note that 
production of neutrinos inside the binary systems has 
been up to now considered in a few different scenarios. 
For example, hadrons accelerated in the jet can interact 
with the radiation of the accretion disk [6|, the stellar 
wind and/or the massive star [H or hadrons accel- 
erated in the pulsar inner magnetosphere or the pulsar 
wind shock can interact with the matter of the massive 
star and/or the accretion disk [Io( ]. More recently, 
production of neutrinos has been considered in the con- 
text of the massive binary systems detected in the TeV 
7-rays, e.g. QUEl]. 



II. PULSAR INSIDE BINARY SYSTEM 

We are the scenario in which the accretion process 
onto the NS surface occurs in the accretor phase. In this 
phase the matter passes through the inner pulsar magne- 
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tosphere onto the surface of the neutron star. It is char- 
acteristic for NSs which are relatively slow rotators with 
a weak surface magnetic field. The gravitational energy 
of the accreting matter is released on the neutron star 
surface, producing a hot spot around the magnetic pole. 
The amount of released energy is re-emitted as thermal 
X-ray emission whose the power can be estimated from: 

L x = GM acc M NS /i? NS « 2 x 10 36 Af 16 ergs" 1 , (1) 

where M acc = 10 16 Afi6 g s _1 is the accretion rate, and 
Lx is the X-ray thermal emission from the polar cap 
region. The radius and the mass of the NS is assumed to 
be i?NS = 10 6 cm and A/ NS = 1.4M . 

The distance from the NS surface at which the mag- 
netic field starts to dominate the dynamics of the in- 
falling matter (the Alfven radius) can be estimated by 
comparing the magnetic field energy density to the ki- 
netic energy density of the wind, B A /8ir — pvf /2, where 
Pa is the magnetic field in the inner neutron star magne- 
tosphere, p = M acc / (irR\vf) is the density of the accret- 
ing matter, V{ — (2GMns/ Pa) 1 / 2 is the free fall velocity 
of the accreting matter, Pa is the Alfven radius, and G 
is the gravitational constant. 

By applying Eq. [TJ and assuming that the magnetic 
field in the neutron star magnetosphere is of the dipole 
type, i.e. Pa = Pns(Pns/Pa) 3 , we estimate the location 
of Ra with respect to the NS surface: 

R A = 7.8 x W 6 Bg /7 M~ 2/7 cm, (2) 

where the magnetic field at the neutron star surface is 
B NS = 10 9 B 9 G. 

Then, we can estimate the magnetic field strength at 
the transition region, 

B A = 3.3 x W 6 Mfl 7 B~ 5/7 G. (3) 

The observed thermal luminosity in X-rays, re-radiated 
from the region of the polar cap, can be calculated from 
Lx = 7ri?p ap oT c 4 ap , where a is the Stefan-Boltzmann con- 
stant. In fact, the emission from the polar cap region is 
well described by a black body spectrum (see e.g. [HI]). 
The radius of the polar cap region on the NS surface, onto 
which the matter falls, and from which thermal X-ray 
radiation is emitted, can be estimated from (assuming a 
dipole structure of the magnetic field) : 

i?ca P = (R 3 NS /R A ) 1/2 « 3.6 x 10 5 £g 2/7 Mil 7 cm. (4) 
Then, the surface temperature of the polar cap has to be, 

T cap = (W^cV) 1/4 « 1-8 x lO^Mjf 8 K. (5) 

In general, the accretion of matter onto the NS can 
occur provided that the radius of the transition region 
(i.e. the Alfven radius Ra) lies inside the light cylinder 
radius of the neutron star, i.e Ra < Rlc = cP/2n, where 
P is the rotational period of the neutron star in seconds 



, and c is the velocity of light. This condition is fulfilled 
for, 

P > P p = 1.7 x l(T 3 Pg /7 M~ 2/7 . (6) 

Moreover, the rotational velocity of the magnetosphere 
at Ra has to be larger than the keplerian velocity of the 
accreting matter. If the rotational velocity, given by 

u rot = 2ttP a /P w 4.8 x 10 7 Bg /7 M w 2/7 /P cm s~\ (7) 
is larger than the keplerian velocity, 

« kop = (GM NS /P A ) 1/2 « 4.8 x W 9 Bs 2/7 Mll 7 cm s -\(8) 

then the matter cannot accrete directly onto the NS sur- 
face. It is partially accumulated close to the transition 
region and partially expelled by the centrifugal force. By 
comparing Eqs. [7] and we get the limiting rotational 
period of the NS below which the matter can accrete onto 
the NS surface: 

P > P a = 0.01B 6 g /7 M w 3/7 . (9) 

NSs with periods within this range, P p and P a , can only 
accrete matter in the propeller scenario. However, for 
periods longer than P a , the accretion process occurs in 
the accretor phase. This last stage is the most interesting 
for this paper. 

III. ACCELERATION OF HADRONS 

As we discussed above, the transition region contains 
turbulent magnetized plasma which provides good con- 
ditions for the acceleration of particles. Both electrons 
and hadrons can be accelerated there. Electrons scat- 
ter on thermal radiation from the polar cap region and 
produce 7-rays by the IC e^- pair cascade process as dis- 
cussed recently in[14|. However, also relativistic hadrons 
can suffer significant energy losses due to radiative pro- 
cesses in collisions with dense thermal radiation from the 
polar cap. These processes produce 7-rays (that undergo 
similar processes as the accelerated electrons) and addi- 
tionally neutrinos that escape from the production re- 
gion without any absorption. The observation of these 
neutrinos by large scale detectors will provide crucial in- 
formation on the relative importance of the acceleration 
of hadrons and electrons in the strongly magnetized and 
turbulent plasma during the accretion process onto the 
NS. 

The maximum power available for the acceleration of 
particles is limited by the extracted energy in the tran- 
sition region. This energy can be supplied by two mech- 
anisms. In the case of a quasi-spherical accretion from 
the stellar wind, the matter has to be accelerated to the 
velocity of the rotating magnetosphere at Pa- The ro- 
tating NS decelerates, providing energy to the turbulent 
region. In the case of accretion through the Lagrangian 
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point, the matter has a large angular momentum, which 
has to be partially lost in the transition region in order 
to guarantee the accretion process up to the NS surface. 
The rotational energy of the accreting matter is then sup- 
plied partially to the transition region and to the NS. As 
a result, the NS reaches the angular momentum and ac- 
celerates. In the first case, the power which has to be 
transfered from the rotating NS to the accreting matter 
can be estimated from: 

£w = Macc« r 2 ot/2 « 10 31 B 8 9 /7 M 3 J 7 P~ 2 ergs" 1 . (10) 

By Using Eq. [9J we can estimate the maximum power 
which can be extracted via accretion from the quasi- 
spherical wind, 

« 10 35 B^ /7 M»l 7 ergs- 1 . (11) 

In the case of accretion through the accretion disk, the 
matter arrives to the transition region with the Keplerian 
velocity. This region is now closer to the NS surface than 
estimated above in Ra (see Eq.[5]) by a factor x ~ 0.1 — 1 
(see (HI). In order to accrete onto the NS surface, the 
matter from the disk has to be slowed down to the rota- 
tional velocity of the NS magnetosphere, i.e. from Wk ep 
to v ro t- Then, the maximum available power extracted 
in the transition region is 

L d « ^cc^Lp - t&t) = ~ ^) ^-(12) 

which gives L d w L™ ax /x for P > P a . 

In conclusion, in both considered here accretion sce- 
narios, from the wind and through the accretion disk, the 
energy is transfered to the turbulent, magnetized plasma 
at Ra- However, in the first case the pulsar loses its rota- 
tional energy but in the second case the pulsar can gain 
rotational energy from the accreting matter and acceler- 
ate its rotation. 

A part, 7], of the power, L w or Ld, can be used for the 
acceleration of hadrons. Below we estimate the charac- 
teristic energies of accelerated hadrons. The acceleration 
rate of hadrons with energy E (and Lorentz factor 7 P ) 
can be parametrized by a simple scaling to the Larmor 
radius of particles within a medium with magnetic field 
B, 

Ace - &E/R h » 1.3 x UfZ_ x M% 7 B; b/7 erg s" 1 ^) 

where £ = lO -1 ^-! is the acceleration parameter, c the 
velocity of light, i?L = E/eB-A the Larmor radius, and e 
electron charge. The maximum energies of the acceler- 
ated hadrons are determined by the balance between the 
acceleration time scale and the time scale which defines 
the confinement of hadrons inside the turbulent region. 
The lower limit on this last time scale gives the escape 
time scale of hadrons from the acceleration site. It can 
be defined as, r osc = R A /v t « 1.1 x 10- 3 Pg /7 M~ 3/T s. 
By comparing this time scale with the acceleration time 



scale, r acc = TO p 7p/P acc (where m p is the proton mass), 
we can estimate the Lorentz factors to which hadrons can 
be accelerated, 

7p « lO^-iB^M^ 7 . (14) 

This estimate does not overcome the maximum possi- 
ble Lorentz factor of particles, which can be accelerated 
within the region with characteristic diameter Ra and 
magnetic field Ba, defined by the condition that the Lar- 
mor radius is smaller than the Alfven radius, i?L < Ra, 
7™ « 6 x 10 6 B 9 1/7 M^ 7 . 

Hadrons with the Lorentz factors estimated in Eq. rjj] 
are above the threshold for pion production in collisions 
with thermal photons from the polar cap region equal to 

7$ » 1.6 x 10 4 P 9 - 1/7 A/r 6 5/28 . 

We can estimate the energy losses of hadrons on col- 
lisions with thermal photons from the NS surface inside 
the transition region, 

P p7 ^ 7r = a P7 cnKE « 1.5 x 10 6 Bg /7 M 33/28 E —,(15) 

where n = n hh (R cllp / ({R A - Rns) + -Reap)) 2 (n « 
nbb-Rcap/Pi f° r -^A 3> -Reap) is the density of the black 
body photons coming from the polar cap at the distance 
of the transition region Ra, and a vl w 3 x lO -28 cm -2 
and K rj 0.15 — 0.35 are the cross section and the in- 
elasticity for pion production due to collision of relativis- 
tic protons with thermal photons. In our estimations we 
apply the average value for K = 0.25. As above, we can 
estimate the maximum possible energies of accelerated 
hadrons due to their energy losses on pion production in 
collisions with thermal photons by comparing Eqs. (13) 
and (15). We obtain the limit on the Lorentz factor of 
hadrons, 

7 P7 _> T « 1.5 x 10 6 £_iB£ /7 Mj£ /28 . (16) 

These maximum Lorentz factors of hadrons are typically 
above the maximum Lorentz factors of hadrons estimated 
above based on their escape from the acceleration region 
especially in the case of large accretion rates (note the 
dependence on the accretion rate in Eqs. (14) and (16). 
Therefore, we conclude that hadrons are at first acceler- 
ated to maximum energies allowed by the escape mecha- 
nism and after that they interact with thermal photons 
during their fall onto the NS surface. Note that, den- 
sity of thermal photons increases towards the NS surface 
with the distance as oc R 2 . Therefore, hadrons interact 
not very far from the acceleration region. 

We consider hadrons which are accelerated with a 
power law or a mono-energetic spectra at the transi- 
tion region. In the first case hadrons with energies be- 
low the threshold for pion production in collisions with 
thermal photons also exist. These hadrons can only 
interact with the matter inside the transition region 
and the matter accreting onto the NS surface. We es- 
timate the density of matter at the transition region, 
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p = Macc/(7ri?|f f ) « 2.2 x 10 16 B^ 6/7 Af^ /r cm" 3 , 
and the interaction rate of hadrons for pion production 
Npp^n = a pp cp(R A /vf) w 2 x 10~ 3 Afi 6 . Therefore, 
hadron-hadron interactions are not efficient at the tran- 
sition region. These lower energy hadrons are convected 
with the in-falling matter on the NS surface. They in- 
teract with much denser matter in the accretion column 
over the polar cap region. However, neutrinos produced 
in hadron-hadron collisions have on average lower ener- 
gies due to large pion multiplicity (see e.g. [16]). We 
show below that these pions (and muons from their de- 
cay) , appearing in a very strong magnetic field on the NS 
surface, can lose energy efficiently via the synchrotron 
process. Therefore, neutrinos from their decay have too 
low energies to produce any observable signal in large 
scale neutrino telescopes. 

IV. PRODUCTION OF NEUTRINOS 

We calculate the spectra of muon neutrinos coming 
from the decay of pions produced in collisions of hadrons 
with the thermal radiation form the NS surface by ap- 
plying a numerical code developed for the interaction of 
hadrons with photons within supernova envelope where 
the condition are quite similar (see [l7| )• Two types of 
hadron injection spectra are considered: (a) a power law 
differential spectrum with the spectral index 2.1 which 
cuts at energies estimated by Eq. [T4j and (b) a mono- 
energetic spectrum with energies given by Eq. 1141 The 
second case is considered as a limiting case of the 
acceleration process of hadrons occurring in the 
reconnection of magnetic field inside the transi- 
tion region. The electric field induced in the re- 
connection region can in principle accelerate par- 
ticles to a very flat spectrum (spectral index lower 
than 2) in which most of the energy is accumu- 
lated with particles at the highest possible ener- 
gies. 

Pions and muons (from their decay) are produced in a 
relatively strong magnetic field at Ra- Therefore, they 
can suffer significant synchrotron energy losses. In or- 
der to check whether or not the synchrotron process can 
change the energies of the produced charged pions and 
muons before their decay, we calculate their synchrotron 
energy loss time scales and compare them with their life 
times. From this comparison, we put an upper limit on 
the Lorentz factor of pions, 

<£ » 7 x l^Bf' 7 M-^\ (17) 

that are able to decay before losing a significant amount 
of energy. We take the effects of synchrotron energy 
losses of pions on the produced spectra of neutrinos in 
the case of their production by hadrons with the Lorentz 
factors above 7® into account by simply replacing their 
7tt by 7*. However, in the case of muons, the critical 
value of the Lorentz factors is about two orders of mag- 
nitude lower. Therefore, neutrinos from muons decays 



TABLE I: Neutrino event rates from a MSP at 5 kpc distance. 



B [G] 


M [g s- 1 ] 


max 


Pa [ms] 


N„ [km 
power 


"V" 1 ] 

mono 


3 x 10 8 


10 17 


2.3 x 10 5 


13 


1.0 


5.1 


10 9 


3 x 10 17 


4.3 x 10 5 


23 


2.1 


5.6 


3 x 10 9 


10 18 


8.4 x 10 5 


36 


5.8 
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have rather low energies. They are not included in our 
calculations of the neutrino event rates. 

For the example calculations, we fix the following pa- 
rameters describing acceleration of hadrons: £_i = 0.1, 
and 77 = 0.1. In table 1, we show the muon neutrino 
event rates that are expected in a km 2 neutrino detector 
from millisecond pulsars at the distance of 5 kpc which 
accrete matter from the stellar wind. The typical pa- 
rameters for the MSPs and the accretion rates are con- 
sidered. The results are shown for a power law spectrum 
of hadrons (marked by power) and a mono-energetic spec- 
trum (marked by mono) . The number of neutrino events 
is estimated by integrating the muon neutrino spectra 
over the probability of their detection, 

where S = 1 km 2 is the surface of the detector, D is the 
distance to the source, P u ^ fJ- (E u ) is the e nerg y dependent 
detection probability of a muon neutrino [201 ] . dN„ / dE„dt 
is the differential neutrino spectrum produced at the 
source, E™ xn — 1 TeV is the minimum energy of pro- 
duced neutrinos over which we integrate the event rate in 
the 1 km 2 neutrino telescope and the maximum energy of 
neutrinos produced in this model is typically in the range 
i?max = 10 4 — 10 5 GeV. This maximum energy of neutri- 
nos can not be simply expressed since it depends on the 
parameters of the model (see Eqs. (14) and (17)). There- 
fore, it is obtained numerically in our calculations. The 
expected neutrino event rates from millisecond pulsars is 
up to a few per km 2 per yr (depending on the model). 
These event rates should be detected significantly by the 
IceCube neutrino detector. 

We also show the neutrino event rates that are ex- 
pected from the two nearby X-ray pulsars inside the bi- 
nary systems GROJ 1744-28 and A 0535+262, which are 
characterized by long rotational periods, strong surface 
magnetic fields, and high accretion rates. Also from these 
sources a few neutrino events per year are expected in an 
IceCube size neutrino detector. 
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TABLE II: Neutrino event rates from X-ray pulsars. 



Name 




L x [erg s 1 ] 


D [KpCJ 


IN ^ [Km yr 










power mono 


GRO 1744-28 


2 x 10 11 


2 x 10 38 


4 


1.6 2.8 


A 0535+262 


4 x 10 12 


9 x 10 37 


1.8 


0.4 0.8 



V. DISCUSSION AND CONCLUSION 

We have considered the model in which hadrons are 
accelerated within the inner pulsar magnetosphere as a 
result of the interaction of accreting matter with the ro- 
tating magnetosphere. The acceleration process occurs 
in the turbulent transition region. Hadrons produce neu- 
trinos mainly in collisions with thermal radiation arriving 
from the hot spot on the NS surface. These hadrons can 
not escape from the NS vicinity since they accrete onto 
the NS surface. Therefore, hadrons accelerated in such 
a model can not contribute to the cosmic ray content 
within the Galaxy. The production process of neutrinos 
is accompanied by the production of 7-ray photons with 
similar spectra and energies. However, these 7-rays ini- 
tiate an Inverse Compton e^- pair cascade in the strong 
thermal radiation from the NS surface (as considered re- 
cently by Bednarek [Oj]). Only 7-rays with MeV-GeV 
energies can escape from the radiation field of the hot 
NS surface. Therefore, we propose that the neutron stars 
accreting at high rates can be observable by the recently 
launched Fermi LAT detector in the GeV 7-rays but also 
by the large scale neutrino detectors such as IceCube and 
KM3NET. 

As an example we estimated the neutrino event rates 
in the km 2 detector for two relatively nearby binary sys- 
tems which accrete at large rate, GRO 1744-2 and A 
0535+262. Other considered sources are expected to 
emit lower neutrino fluxes due to the larger distances 



or lower accretion rates. We also estimate the neutrino 
event rates expected from the millisecond pulsars at the 
typical distance of globular clusters. It is clear (from 
Table 1) that some accreting millisecond pulsars could 
become detectable neutrino sources. 

Different factors can enhance or reduce the neutrino 
event rates reported in Tables 1 and 2. The obvious ones 
are related to the distance to the source and the variable 
neutrino emission due to a change in the accretion rate 
caused e.g. by the eccentric orbit of the neutron star. 
Moreover, in the case of the accretion process through the 
accretion disk, the power available for the acceleration 
of particles can be enhanced by a factor of x _1 (\ is 
typically in the range ~ 0.1 — 1, see Eq. 12) with respect 
to the case of the quasi-spherical accretion. Therefore, 
the neutrino event rates reported in Table 1 should be 
scaled as well by that factor. From another site, most 
of the observed millisecond pulsars belong to globular- 
clusters. It is expected that a massive globular cluster 
can contain up to ~ 100 MSPs[21]. In fact, inside the 
globular clusters Tuc 47 and Ter 5, already ~ 20-30 MSP 
have been discovered A significant amount of the 
MSP cannot be discovered directly by their pulsed radio 
emission since they are inside compact binary systems 
(so called hidden MSPsf22]). These hidden MSP can 
become neutrino sources via the production mechanisms, 
that have being discussed in our model. Therefore, the 
neutrino event rates expected from the whole population 
of MSPs in a specific globular cluster can be enhanced 
due to the cumulative contribution from many sources. 
On the other hand, non-observation of neutrinos from 
the globular clusters give important constraint on the 
number of hidden MSPs in these objects. 
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